Correlation in superradiance: A closed-form approach to cooperative effects 
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We have developed a novel method to describe superradiance and related cooperative and collec- 
tive effects in a closed form. Using the method we derive a two-atom master equation in which any 
complexity of atomic levels, semiclassical coupling fields and quantum fluctuations in the fields can 
be included, at least in principle. As an example, we consider the dynamics of an initially inverted 
two-level system and show how even such in a simple system phenomena such as the initial radiation 
burst or broadening due to dipole-dipole interactions occur, but it is also possible to estimate the 
population of the subradiant state during the radiative decay. Finally, we find that correlation only, 
not entanglement is responsible for superradiance. 



I. INTRODUCTION 

Following the appearance of Dicke's paper [| on spon- 
taneous emission of radiation by an assembly of atoms, 
cooperative effects such as superradiance were studied by 
many authors, e.g. Q. Recently, there was a revived in- 
terest in, so called, Dicke states because of their potential 
application in quantum information processing 0- [3- 13- EJ 
and their importance in the behavior of Bose-Einstein 
Condensates 0. 

The main reasons to revisit superradiance in this arti- 
cle are (a) the need for an improved calculational model 
to support the revived interest, and (b) a set of novel phe- 
nomena that could be discovered using this model, such 
as subradiance and chirping, (c) the need for a simpli- 
fied formalism to treat more complicated level systems, 
and (d) the improved understanding of the phenomena 
of collective vs. cooperative effects. 

We have recently developed a novel method to study 
collective effects in optically dense media @. Whereas 
collective effects arc caused by a radiative dipole-dipole 
interaction between the atoms, cooperative effects also 
include the exchange of virtual photons between the 
atoms that leads to the formation of Dicke states. The 
difficulty of calculating effects including atom-atom co- 
operation has been the intractably large number of in- 
terconnected degrees of freedom, even if just a few par- 
ticles are involved. In this article we incorporate coop- 
erative effects into our formalism, leading to an effective 
two-atom master equation formalism for superradiance. 
From the simplest version of this calculation we are able 
to obtain good agreement with observed signatures of 
superradiance including the effects of dissipation and the 
unique temporal build-up of a sharp flash of radiation. 
In addition, our new formalism allows more complicated 
level structures, additional fields, and polarization to be 
easily considered as well. 

What is the main difference between collective and 
cooperative effects? In this context, we define as "col- 
lective effects" all those, where a multi-particle ensem- 
ble responds differently than just the sum of all single- 



particle responses. Examples of collective effects in a 
atomic gas with radiative interactions include amplified 
spontaneous emission and radiation trapping (see e.g. 
[8| )■ "Cooperative effects," on the other hand, specif- 
ically arise from the exchange interaction between any 
two particles. The best-known example for this exchange 
interaction is the energy splitting between the symmet- 
ric and anti-symmetric two-atom superposition states, 
(\ab) ± |6a))/-\/2, where |a) and \b) are the excited and 
ground states of a single two-level atom. In this arti- 
cle, we will quantify the exchange interaction as the non- 
diagonal two-atom term of a two-atom density matrix 
p, Trp|afe) (ba\ = p a b,ba- We will show that the transi- 
tion between non-cooperative and cooperative collective 
effects is given by the transition between p a b,ba = Q tp 
Pab.ba 7^ 0. All typical signatures of superradiant behav- 
ior can be explained using this term. 
Since the term Pab.ba 

describes the coherence between 
any two atoms, we will refer to it as "correlation," and 
we will show (see, e.g., Fig-EJ) that the correlation is the 
reason for the build-up of superradiance. 



II. MODEL 

We consider an arrangement of N atoms which mu- 
tually interact via quantized electric field in dipol and 
rotating-wave approximations (similar as was done in 
0)- We are interested in cooperative phenomena of the 
ensemble so we distinguish two atoms in the arrange- 
ment, say 1 and 2, and write the total Hamiltonian as 

H = H md +J2{ H o- & + £i ) } 

i=l,2 

+ E { H i-Pi(z j+ £ )), (i) 

for the probe atoms i = 1, 2 and all other atoms j ^ 1,2. 
Here, Hq J and Hfi e id are the free Hamiltonian of the i,j- 
th atom at location ?i ,j , and that of the quantized radi- 
ation field, E, respectively. Furthermore, Ei.j = E(fij) 



2 



is the quantized field at the position of the i, j-th atom, 
while £j = £{rj) is the external classical driving field. 
Pi j — p(fij) is the atom dipole operator. 



for the /J, = x, y, z-components along each of the coor- 
dinate axes. Thus, with a rotating wave approximation, 
Eq. J2J becomes 



3- 



sf 



T c exp 



{(iw) + (i( 2 ))} , 



(4) 



FIG. 1: The integration on the Schwinger-Keldysh contour 
follows the dark line from — oo via t to — oo. The time ordering 
is normal time ordering on the upper branch, inverse time 
ordering on the lower branch, and times on the lower come 
before times on the upper branch. 

In the next step we move to an interaction picture us- 
ing the time evolution operator on the appropriately cho- 
sen Schwinger-Keldysh contour shown in Fig.lTTl 



Sc = T c exp 




(2) 



where Tc is a time ordering operator on the Schwinger- 
Keldysh contour, and is identical to standard time or- 
dering operator T on the upper or + branch, and to its 
inverse T _1 on the lower or — branch. By this formal 
procedure each physical time t is replaced by two times 
on the contour, t + when going from — oo to +oo on the 
upper or + branch, and i_ when going in opposite di- 
rection on the lower or — branch. In this notation, an 
infinitesimal complex quantity is added to what was pre- 
viously real time r in order to obtain f, where the sign 
of the infinitesimal complex part depends on the branch 
of the Schwinger-Keldysh contour. 

We find the expectation value of Eq. @ by replacing 
formally 



T c exp{i}) = cxpj]r((T c i" 



(3) 



where A is the operator in the exponent of Eq. @ . The 
cumulants contain the pure correlation, e.g., ((x)} = (x), 
((xy)) = (xy) — (x) (y), etc. Assuming that the radiation 
field is Gaussian allows us to bring to a tractable form 
the cumulants involving the radiation field, which appear 
when averaging the evolution operator [j| 0, ^| with 
respect to the radiation field. This told, in the effective 
evolution operator only terms up to quadratic in electric 
field survive. 

In addition,we separate positive and negative fre- 
quency components of the atom dipole operator p M and 



the electric fields and £ M , 



where 

id)' 

1(2) 



c i=1 > 2 

7j I dfi / df 2 > 

2h 2 J J f? 12 

C C »>J— M 



where we used £l = £ + \Ej for the local field seen by 

each of the two atoms. The two cumulants that appear 
in Eq. (QJ are 

A„J* (n , r 2 ) = {{T C E^ (n)E+ (f 2 )» e-^-^a) 
CW(ri,f 2 ) = ((TcEf^Ej^)))^-^ )(5b) 

The function D, for example, can then be expressed as 



AAt,j>(n,f 2 ) = c 



'CJi' ' 



(6) 



with 



D 



//.;/ 



((TE-(r f i ,T 1+ )E+(? j ,T 2+ ))) 
E-(r l ,T l _)E+(r ]1 T 2+ ))^ = D~+ , 



'E+(r ] ,T2-)E-(r l ,r 1+ ))] 



[IV ' 



where we used the properties of the time ordering oper- 
ator Tc on the Schwinger-Keldysh contour. In the ex- 
pression above a subscript sign "+" or "-" to the time 
variable ti,2 indicates the branch of the Keldysh contour 
on which the time variable is located. Similar expression 
can be written down for function C as well. We assume 
that there is only one type of atoms so we set 



Mr) 



E+(f) + E-(f), 



Pu(f)e il 



p£{n,t) 



(7) 



E+(f)e-^ + E-(f)^, 



where p is the dipole matrix element for all polarizations, 
_ _.^ t and <?i(£| is the dimensionless dipole vector operator of 

= £(i (t) + £p \ T )-i^n {^) e + Sjidfivldual' atom i. For real times, the effective time evo- 
lution operator Q becomes, after all the contributions 
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from the two branches are added together, 
t +00 



c - T c exp \ \ dr 



Sf 



E I E KfJ v( T «) f L, v ( T «) + H -' 



i=l,2 re= 
+00 +00 



2?i 2 



C?Ti C?T2 



(8) 



i,j'=l,2 re,A=+,- 
«A 



are real and imaginery part of the Fourier transformed 
cumulant D '", and so they are mutually related via the 
Kramers-Kronig relationship. 

Term T describes the decay and pump rates induced by 
the incoherent photons inside the medium, while the term 
7 is the spontaneous "down rate" in the atomic medium. 
Term H describes a collective light shift in addition to 
the inhomogeneous broadening, and is incorporated in 
our formalism. Term h is a spontaneous contribution to 
the light shift, the diagonal terms of which amount to 
the Lamb shift. As we are not interested in the Lamb 
shift, we consider these terms either included in the free 
Hamiltonian Hq (diagonal), or zero (off-diagonal). With 
the help of Eq. 11UII we can write the effective two-atom 
density matrix equation as 



Here, the superscripts in radiation field averages C and 
D determine on which of the branches of the Schwinger- 
Keldysh contour respective time argument lies, example 
of which was given in Eq. @ . Additionally, we assume 
that the field coherence is shorter in duration than the 
coherence of the atomic operators cr,a\ e.g., 

J dTidT 2 u{Ti) ■ F(ri,r 2 ) • cr t (r 2 ) 

= J dr J dT'a(T- T -)F{T- T -,T+ T -)a\r+ T -) 

« J dTo{T){J dT'F{T- T -,T+ T -)YHr) (9) 

(Itcf(t)!F(t, uj)a^ (r). 



This yields a more convenient way to reexpress the quan- 
tum correction functions C kX and D kX , with k, A £ 
{+,— I of Eq. (J8J, based on the following functions (see 
also for details), 



+00 



P^Pv 



dr' 



(10a) 



EVP" I dT > 



fl f /Li,, jv(T,Ul') 

Hin dv {T,u>) = — V \ dJ lfl ^ K ' ; 
p J 2tt J u - u' 

—00 

+00 

hi^,jv{T,uj) = —V / duj 



(10b) 

<[^(r),^ J v(r + T')]>e^' > 
(10c) 



UI — to' 



(lOd) 



Here, V denotes the principal value of the integral that 
follows it. Please note T, Eq. (|10a|l . and H, Eq. (|10c|) 



[Ho,p] 



(11) 



j=l,2 



H,v .3=1,2 



r 

x 11x.1v 



E E 

i,j=l,2 
i,j=l,2 



where £l,^(^) is the local (Lorentz-Lorenz) field felt by 
the j-th atom. Eq. I|ll|l is one of the main results of 
this paper. Up to this point we assumed a homogenous 
atomic medium, made a Markov approximation, and ne- 
glected correlations of higher than second order. All co- 
herent fields are displayed in the term together with 
the first-order quantum corrections which take care of 
the local field. All second-order corrections are in the 
terms with H and L. This equation looks similar to 
the effective master equation for large times, derived in 
Refs. [HHii. Note, however, that Eq. 

is a two-atom 

equation! In addition, the second-order quantum correc- 
tion terms, which have quite a different form in our case, 
contain most of the physics in the context of cooperative 
effects. It is also possible to find a closed-form expression 
for them, as will be seen in the next section. 



III. COOPERATIVE EFFECTS IN A 
HOMOGENEOUS GAS OF TWO-LEVEL ATOMS 



For a gas of two-level atoms there is only one transi- 
tion, and the coordinate indices /x and v can be dropped; 
the effective two-atom two-level master equation <|11[) as- 



4 



sumes the form 
i 



[Ho,p]+ J2 ip cr J f i(^)+ '] f i( f S)>P 



H 



h 

1 y 

3=1,2 
i,3'=l,2 

1 



3 = 1,2 



3,3 [[ Cr 3' CT ] 



(12) 



~0 H ( T i,3+7i,j) 



<7j , <?iP 



i,3'=l,2 



Let us introduce the notation for the two-atom density 
matrix as following, 



Pa(3,~t& = (ai72|p|/M2 



(13) 



Each of the two atoms can be in the ground state \b), or 
the excited state |o), and 1 0:1,72) is a product \a) 1 ■ |7) 2 
of the states of the first and second atom. The atomic 
lowering operator a assumes the form a = \b) (a\. This 
allows us to introduce three real functions, namely a, d 
and n, as follows 



2 a — 2p aa ^ aa + p aa ,bb + Pbb,ao 
d = Paa,bb — Pbb,aa, 

Paa,aa Paa,bb Pbb,aa 



(14a) 
(14b) 
Pbb,bb- (14c) 



Thus, a{f\,r2]t) is the single-atom excited-state popula- 
tion averaged over both atoms, d{f\,r2",t) is the differ- 
ence in excited-state population for the two atoms, and 
n{fi,?2',t) is the product of the inversions of the two 
atoms. In new variables the equations of motion read 



a = -(2r + + 7 )a + r + - T_d , 
d = -(2r + + 7 )d-2r_(2a- 1) 
h = -2(2r+ + 7)n-2 7 (2a-l) 
+4 (f + +f_)p ah . ba +4(f + 

Pab,ba = - (7 + 2r + + 1A12) Pab.ba 

+ (f+ + f_)n. 
where we used 

2r± = r n ±r 22 

2f± = r 12 ±r 21 
2 

Al 2 = — (#11 — H 2 2)- 

n 



-r_ 



(15a) 
(15b) 

' Pba.kl^c) 

(15d) 
(15e) 



(16a) 
(16b) 

(16c) 



It will become clear, that the cross-coupling or correla- 
tion term p a b,ba carries the cooperative physics. 



nate dependence can be dropped from the variables de- 
scribing the system, a, n, d, and p a b,ab- That is, we ne- 
glect retardation effects of propagation of the electromag- 
netic field through the sample, that is, all time changes in 
the atomic variables propagate instantaneously through 
the sample. This is a good approximation as long as the 
"cooperative" time tq = (A/" CX c7m) _1 , with A/" eX c the den- 
sity of atoms in the excited state, 7 the vacuum decay, 
and p a geometric factor, is longer than the maximum 
propagation time • 

In the small sample approximation, thus, there is no 
difference between the atoms 1 and 2 with respect to 
their decay rates. That is, in Eq. I|15|) we set T + = Tn = 
r 22 = r. Similarly, there is no spatial difference between 
the decay rates so T_ = f _ = 0. For other decay rates 
we have T + = Ti 2 = T 2 i = T. Finally, because in a 
homogeneous sample the result is invariant with respect 
to exchange of the indices of two atoms, we have Ai 2 = 0. 
This said, the system (fT5|) reduces to, 



a 
h 

Pab,ba 



-(2r + 7 ) a + i 

-2(2r + 7 )n- 

-(2r + j)pab,bc 



27(2a - 
+ fn. 



(17a) 
8f p aMQ (17b) 
(17c) 



Let us now discuss the decay matrices Tij, where 
i, j = 1,2, and which we introduced in Eq. (|10fl . in greater 
detail. Our discussion is based on the formalism for 
calculation the quantum corrections to the electric field 
in a responsive media, that we have previously devel- 
oped in [8|. We need the exact Green's function for the 
Maxwelf field E, 



D^(i,2) = ((T c E-(i)Et(2))) 



(18) 



where Eq would be the field without medium present, as 
calculated along the Schwinger-Keldysh contour. Here, 
we abbreviate ?x,ii by 1, etc. The Green's function for 
the free field Eq is given by 



D^ V {1,2) = {{T c E^(l)E+ v {2))) 



(19) 



The contour Green's function D contains four real-time 
functions, depending where on the Keldysh-Schwinger 
contour the two time points are located, D ++ , D h , 
, and D , as discussed in 14]. We find approxi- 
mately |3, D+ + w D Mv , D, 
and Dq w —D 
and advanced Green's functions, respectively. A solution 
to the atom-field interaction can be written in terms of 
a Dyson equation, 



~0,D+- 
where D ret and i5 adv are retarded 



rjadv 
u 



n rct 



A. Quantum corrections in a small sample 
approximation 



II dt ' ldt ' 2 



(20) 



We further discuss the Eq. I|15|l in terms of small sam- 
ple approximation, where we assume that all the coordi- 



-.3 . -.3 . 



V1V2 



d?ldP" 2 D%(l, 2')^ v (2', 2), 
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where 



+00 

D%{1,2)=DI%{1,2)- JJ dt' x dt 



' "'" l 2 



(21) 



df^D^l, l')P^(l', 2')Dtf(2', 2), 

V1V2 

with Z)^ v (l,2) = Z>02,1). In Eqs. fZQ and fig- 
ure two source (polarization) functions. In the lowest 
order in the atom-field coupling, the polarization func- 
tion is given by a correlation function of dipole operators 
of noninteracting atoms, 



(22) 

and for retarded polarization, 

P^ c \r u r j] t l M)= (23) 

In Eqs. i|22|) and 1|23|) the superscript 1 stands for a one- 
atom source function, i.e., i — j, while 2 is for a two-atom 
source function, and i 7^ j. The cumulants can be found 
by using the quantum-regression theorem and a Laplace 
transformation, as done below. 

The differential equation for the single atom density- 
matrix cross-term, p a b, is found by tracing over the sec- 
ond atom, 

Pab = Pab,aa + PabM = ~ ( ^ + T + iA^j p ab (24) 

with A = Aj = 2Hu/h. Upon comparison with the same 
equation for a dilute gas, we find that the second-order 
correction Ha in this case takes on the role of a detuning 
from resonance. Since, in general, Ha, along with T, 
changes over time, a chirp of the radiated light can be 
expected. 

We can Laplace-transform the source functions from 
r to A with t\ — > t + t and t 2 —* t for t\ > t 2 (and 
vice versa). Thus we find for Eq. (|24|l (and its complex 
conjugate) 



Pab = 



Pab{r = 0) 
A + 7 /2 + r + iA : 



Pba = 



Pba(T = 0) 



A + 7/2 + T-lA' 

(25) 

Here (crj(A, i)) = p a b- Using these solutions and the 
quantum-regression theorem |l5j , we find for the Laplace- 
transformed single-atom cumulants (i = j) 



oft) 



2+r-iA 
1— aft) 



A+7/2+r+iA 
and for the two-atom cumulants (i =/= j) 
aJ(A,^(t))) = x+ Jt ] r _ lA = 
^(A,t)<xt(i)\\ = A+7 / 2 ( + ) r+lA = 



(26) 



4(*)cri(A,t) 
a 3 {t)a\{\t) 



(27) 



In a simple case like this we can go from Laplace- to 
Fourier-transform by just replacing A by iA', where A' 
is a frequency. Then, we find for the source functions in 
Fourier space 



P^ s (x,f;A,t) 
P (2)s (f,r; A,t) 
P^ Ict (x,f;A,t) 



2a(x, t)(7/2 + T) 



h 2 



®-M 2 

h 2 (7/2 + r) 2 + (A' - A) 



(7/2 + r) 2 + (A' - A) 2 
2x(x, t)( 7 /2 + r) 



<28) 



1 -2a(x,t) 



h 2 7/2 + r + i(A' - A) ' 



(30) 



where fi — fj = x and fj — f. Since in this case there is 
no external light field present, we can set A' = 0. 

The equations above allow us to write down closed 
form expressions for both T and F, where 



r(A,f) 



n Jv 



Ti 

dx 



f(A,i) 



2 D-+{n =f 2 ,A,t) = (31a) 
L> rct (r,£,A,i) 2 p( ls \x,A 7 t) 
dx 1 x 2 b Tet {f,x 1 ,A,t) 

2,A,t) 

(31b) 



h j jy 

D* Iet (r,x 2 ,A,t) P^(x u x 2 ,A,t) 
o 2 ~ 

^D-+(f 1 ^r 2 ,A,t) = 



h — \ dxD ret (fi,f,A,i) 

h Jv 



D" et (f 2 ,x,A,t) P^(x,A,t) + 



dx 1 x 2 D Tet {f 1 ,x 1 , A,t) 



n j Jv 

D* Ict (f 2 ,x 2 ,A,t)P^(x 1 ,x 2 ,A,t) 



where we use the Fourier transforms for all the functions. 
The A' now denotes the detuning from resonance of the 
particular Fourier components at frequency lj — lj^+A', 
where uj a b is the atomic resonance frequency. In or- 
der to solve Eqs. I|31|) we make an additional approxima- 
tion where we neglect the coordinate dependence of the 
atomic variables a,n, p a b,ba on fx,f% in the integration. 
This is justified if we assume a much weaker coordinate 
dependence of the atomic dynamics than of the field cor- 
relations. In Ref. |8| we have discussed in great detail a 
way of calculating the kernel D ret {x\ , x 2 , A,t). Here we 
just give the final result 



D^{x u x 2 .A:t) 



-i go (A) r 



(32) 



6 7T £o C 2 

where r = \xi — x%\. Here, an average over all dipole 
directions is taken. The wave vector go is given by 

Tim 



<?o(A) = - 1 



c \ 3 £o c 



P^ Tet {r,A;t) 



(33) 



In Eqs. IMl'li and l|33|) . ui stands for the frequency of the 
emitted light. 
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We then obtain a self-consistent expression for the 
quantum corrections as follows 



a(t) 



2a(t) - 1 



+2 7 CV 



2C(A) _ 



(34a) 



Y) + A* 



PabM (t)I(C(A),g(A)) 



r = 3jCg 



r} + A* 



a(t)I(C(A),g(A)) 



(34b) 



-2 7 CV lTf 



Tj + A* 



P„i )6 a(t)/(C(A),fif(A)) 



In the expressions above both decay rates, T and F, them- 
selves implicitely depend on the atomic variables a,n,d, 
and p a b,ba of Eq. I|14|) . so they have to be self-consistently 
calculated. Different quantities that appear in Eq. i|34|) 
are given as follows, 



r 



/ 



= 2 + r 



C(A) 



Co = fc e Z(2a(i)-l), g(A) = ?- 2 |c(A) 

Q = 



_ _d 

A 2c ' 2 2 

((C— i)c c +cos e) +(^0^— sin g) 
(C 2 +e 2 )" 



J(C,<?) = 

(35) 

We see that the Eqs. (|34|l and (|35|l rely on two parame- 
ters: effective density, which is also known as the coop- 
erativity parameter, C, given by 



C = 



2ttc 3 M 



(36) 



which is proportional to the total number of atoms in a 
cubic wavelength of the emitted radiation, where AT is 
the density of the atoms; and the effective radial size, g, 
Eq. H35fl . proportional to the ratio of diameter d of the 
sample to the wave length A of the emitted radiation. 

Now, the role of A = 2Hu/fr becomes clear: As a real 
part of the dispersion, A is related to imaginery part T, 
a decay rate due to cooperativity effects, via a Kramers- 
Kronig relationship 



A = —V 

IT 



dA 



, T(a,x,T,A') 



(37) 



that in our modeling we have neglected the time retar- 
dation effects within the sample. When detecting the 
emitted radiation one would have to add some delay be- 
fore the radition reaches the detectors. 

We now examine short time behavior of the system 
(|17|l when all the atoms are initially in the excited state, 
\a). The values for the atomic variables a(0), x(0) and 
n(0) are then a(0) = l,n(0) = 1, p a b,ba — 0. While the 
system of equations (fTTjl provides a way of calculating the 
chirp A from the spectral distribution V = T(a, n, x, A), 
for simplicity we set A = for rest of the discussed cal- 
culations. We choose as values for the effective density 
C = 10 and the effective size g = 10. Different aspects 
of the solution of system l|17|) . given the initial values 
and relevant parameters as described above, is shown in 
Fig-El (Top panel) reveals a typical signature of the co- 
operative enhancement to the single atom emision rate at 
the early stages of deexcitation. (Middle panel) depicts 
the quantities p a b,ba = Pab,ba( T ) anc ^ n = n ( T )- Here, the 
most striking is the behavior of the two-atom entangle- 
ment, p a b,ba, which grows from zero to some maximum 
value before dropping back to zero (for times <C 7 , not 
shown here). This feature, a build-up of entanglement 
with a significant increase in emitted intensity per atom 
compared to a single atom, is a characterteristic signa- 
ture of superradiance, as discussed in, e.g., Ref. [T(| . 

One of the features of the superradiance is that the 
intensity per atom depends strongly on the density of an 
optical medium under consideration. In our theory the 
measure of density is the effective density C, Eq. Q36p. 
We perform a series of calculations where we vary C, and 
calculate max{— a(r)} r = — d(r max ) and r max for each of 
the solutions. The results of this calculation are shown in 
Tabled maxima of the decay rate are proportional to the 
effective densityC, while the times at which maxima occur 
depend are reversely proportional to the effective density 
C. This is in good agreement with the expected super- 
radiance scaling, that is, the emission per atom — <z(t) is 
indeed proportional to C, while the times at which the 
maxima are obtained scale like 1/C. The curves with the 



c 


max{— d} T 


Tmax (in 7 L ) 


10 


57 


0.0018 


20 


115 


0.0008 


30 


172 


0.0006 



In a two-level system, A is very small, and we can set 
it self-consistently to zero in what follows. The detailed 
format of A will be discussed in an upcoming publication. 



B. Cooperative Phenomena 

A physical measure of the dynamics in the system is the 
average intensity of the emitted radiation, rather than 
the effective population of the excited state. The in- 
tensity of radiation is directly proportional to the term 
— a(r) by energy conservation. One has to keep in mind 



TABLE I: Dependence of the maximum of the intensity per 
atom — a and the scaled time r max , at which the maximum 
occured, as a function of C. The quantity — d is a linear 
function of C, and thus confirming the superradiant nature 
of the decay. The times r ma x at which the maximum of the 
radiation is achieved, scale like 1/C. 

calculated intensity per single atom for the cooperativity 
parameter C = 10, 20, 30 are given in Fig.|3J 

In the original Dicke's paper 0,the superradiance oc- 
curs through a decay over the subset of symmetric many- 
atom states. In our effective two-atom two-level model, 
we may introduce two alternative matrix elements: a 
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1 0.8 

% 0.6 

I 0.4 
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0.2 








a(T) 

. . . n(l) 












0.005 0.01 0.015 

Scaled Time x ( yt ) 



0.02 



FIG. 2: Cooperative effects in opticaliy dense medium at 



short time scaies (< 7 



(Top panei) Intensity per 



atom, — d(r), of spontaneous radiation is at the initial mo- 
ments much stronger than that of a single atom. (Middle 
panel) Behavior of the atomic variables shows an expected 
decay of the excited state, a(r), and the population differ- 
ence, n(r)), with a slow build-up of entanglement, p a b,ba(j), 
from initial zero value. (Lower panel) Self-consistently calcu- 
lated F and r along the trajectory of the system (1171 in phase 
space. Found values are two orders of magnitude larger than 
the spontaneous decay rate 7 (here 7=1). 




0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 0.009 0.01 
Scaled Time Tiny! 

FIG. 3: Superradiant intensity per atom, — d(r), in dense op- 
tical media of cooperativity parameter C = 10, 20, 30. Greater 
the cooperativity parameter, greater the maximum momen- 
tary deexcitation rate T = F(r), see also TableQ 



subradiant, p , for the population of anti-symmetric 

two- atom state, and a superradiant, for the sym- 

metric two- atom state, where |±) = \/\/2{\ab) ± \ba}). 
This is permissible if the size of the atomic sample is 
small, because the effects of change of phase due to 
propagation can be neglected. In terms of atomic vari- 
ables a(r),n(r) and p a b,ba( T )' these matrix elements are 
p±± = 1/4(1— n(r)±4 p a i Jj f, (T)), The behavior of the ma- 
trix elements under the conditions exhibiting superradi- 
ant behavior, Eq. I|17|) . and with the cooperativity param- 
eter C = 10 is shown in Fig.0] Here we observe that the 
dipole-dipole interaction, which is built into our model, 
indeed contributes to populating the anti-symmetric two- 
atom state, whereas the population of |+) and | — ) equal- 
izes for long times when decoherence overrules correla- 
tion. This feature is absent in standard treatments of 
superradiation, e.g. |l6j . 




0.0001 0.001 



0.01 0.1 
Scaled Time T in y t 



FIG. 4: Superradiant, p++, and subradiant p — density ma- 
trix element in logarithmic time for a dense optical media 
with C = 10. The solution of system (I17H shows that in a 
superradiant decay both, superradiant and subradiant, states 
get populated. As expected, there is a significant difference 
between the two, as long as there is an entanglement present, 

pabM^) / 0. 



IV. CORRELATION AND ENTANGLEMENT 

In this paragraph we would like to clarify the respec- 
tive roles of crrclation and entanglement in superradi- 
ance. Correlation, i.e., a two-atom cooperative effect, 
is stored in the non-diagonal two- atom matrix element, 
PabMi an d can be measured as the energy splitting be- 
tween the |+) and |— ) state. This quantity clearly is 
responsible for the build-up of superradiation (see, e.g., 
Eqs. {13), {T7J|, and (EE})- Entanglement, however, de- 
fined to be a quantity describing how far the system is 
away from a product state p pro d — J2i Pi^ ® Pi i where 
Pi^ are the single atom density matrices. Wootters 
has found a measure for the entanglement of two systems. 
It turns out that in all two- atom superradiant systems, 
even those with a maximum grade of coherence, there is 
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no entanglement present. The proof is easy, since p a bba 
always has to be smaller or equal to ^p aa pbb, where p aa 
and pbb are the single-atom populations. 

Obviously, at this time it is not possible to make any 
final statement about the entanglement of more than two 
atoms with cooperative interactions, but since we treat 
only effects up to second order in the atomic interaction 
exactly, it is quite clear that also in our system there is 
no entanglement present. In conclusion, supcrradiance is 
an effect based on correlation, but not on entanglement. 

V. CONCLUSION 

We have presented a method for modeling the coop- 
erative effects in a dense optical medium, based on an 
effective two-atom two-level model. In its simplest ver- 
sion, a small sample approximation, we are able to obtain 
a system of equations which describes the evolution of a 



two-atom two-level density matrix. Solving the system 
for an initially inverted system, we are able to predict 
the cooperative behavior with all the features of super- 
radiance: coherent build-up of entanglement, with the 
decay rate up to two orders of magnitude greater than 
the spontaneous decay rate. Additionally, the maximum 
of the decay rate T scales linearly with the number of 
particles indicating that the maximum intensity of the 
deexcitation pulse scales as N 2 , the square of the num- 
ber of particles. Our model requires two parameters, the 
cooperativity parameter C, and the effective size of the 
system in terms of the emitted wavelength. Theory pre- 
dicts two novel phenomena: chirp, or the change in tran- 
sition frequency due to the cooperative phenomena and 
subradiance, both of which are outside a standard de- 
scription of superradiance. Accounting for the shape of 
the optical medium in the formalism, and in particular 
analysis of cigar-shaped atomic sample, will be subject 
of future work. 
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